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We consider particle creation (the Dynamical Casimir effect) in a uniformly contracting ideal 
one-dimensional cavity non-perturbatively. The exact expression for the energy spectrum of created 
particles is obtained and its dependence on parameters of the problem is discussed. Unexpectedly, 
the number of created particles depends on the duration of the cavity contracting non-monotonously. 
This is explained by quantum interference of the events of particle creation which are taking place 
only at the moments of acceleration and deceleration of a boundary, while stable particle states exist 
(and thus no particles are created) at the time of contracting. 

PACS numbers: 42.50.Ct,42.50.Dv,03.65.-w 



I. INTRODUCTION 



j^ ' During the last two decades the dynamical Casimir effect (DCE), the effect of photon creation in an empty nonsta- 

^ tionary cavity [l|, attracted considerable attention in literature. Being a firm prediction of quantum field theory, the 
O"^' DCE has not been observed yet experimentally. This is probably the main reason for such interest to it, and there- 
J> I fore search for a realistic scheme for experimental observation of the DCE is an actual problem of modern quantum 
•"^ ■ electrodynamics. 

K% ' A nonstationary cavity can be realized by two possible ways: (i) a cavity with moving boundaries, or (ii) a cavity 

t^ with fixed shape but with varying boundary conditions. In case (i) a significant number of photons could be created 

■ ■ ■ ' if the velocity of boundaries is close to the speed of light. The corresponding experimental realization would be 

obviously a difficult technical problem. Alternatively, it was suggested to utilize the vibrating cavity and accumulate 

the effect by tuning the frequency of small mechanical oscillations of the boundaries (which could be inducedjC.g., 

by a sound wave [3, or due to piezoelectric effect Q) in resonance with an eigenfrequency of the cavity [1, 1MB 0131' 

see also the recent review j^- In this realization, an extremely high accuracy of frequency tuning must be provided. 

Owing to amazing progress in modern laser technology and semiconductor electronics, the variant (ii) could be 

considered as a more realistic one [13, [HI [l3 ■ The optical properties of a semiconductor film located on a dielectric 

base can be changed, e.g., by means of electron-hole plasma creation by a strong femtosecond laser pulse [iJl, or due 

to injection of carriers by a powerful electric pulse [I3j . In both cases, variation of optical properties of the boundary 

can be fast enough and well controlled. 

However, theoretical investigation of the first variant is simpler from mathematical point of view, at least in the 
framework of a one-dimensional model. Sometimes in the framework of the first approach exact solutions to the 
problem can be found for some special laws of motion. Besides, one may hope that very fast displacement of a 
boundary acts like " almost instant" creation of a new boundary due to changing of optical properties of the medium 
inside the cavity. Therefore we will consider a cavity with moving boundaries in this paper, and will be especially 
interested in the case of ultrarelativistic motion. 

We will present a new exact solution for the DCE problem in a one-dimensional ideal cavity uniformly contracting 



during the time interval < t < T and otherwise stationary. A complete set of solutions for a classical field inside a 
one- or even three-dimensional cavity at the time of uniform contracting was known for years, see Ref. |2|- It was used, 
for example, in Ref. [liL Il5l|. for calculation of the Casimir force acting between the boundaries of a relativistically 
squeezing or expanding cavity. In contrast, our solution for a one-dimensional quantum field in a uniformly contracting 
cavity seems to be unknown in literature, and allows one to study the DCE for this particular case non-perturbativly. 
We have discovered a new and unexpected effect of periodic dependence of the number of created particles on some 
variable Xf determined by the time of contracting T and velocity of the boundary v. It is explained by quantum 
interference of the events of particle creation which are taking place only at the moments of acceleration (t = 0) and 
deceleration {t = T) of a boundary, while stable particle states exist (and thus no particles are created) at the time 
of contracting. 

We review the Hamiltonian approach in Section m Consideration of a field quantization procedure in a uniformly 
contracting cavity is presented in Section UTTI Section Hvl is devoted to calculation of the energy spectrum of particles 
created in the cavity. The periodic dependence of all measurable quantities on Xf is indicated, and the optimal 
conditions for the DCE are derived. The discussion of the results and conclusions are given in Section We use 
natural units h = c = I throughout the paper. 

II. HAMILTONIAN AND OTHER PRELIMINARIES 

We will use the Hamiltonian approach which in application to DCE was apparently first formulated in Ref. |l6| . 
see also l(|, |9| . In this section we begin with derivation of Hamiltonian for a one-dimensional problem based on the 
well known properties of time-dependent canonical transformations in classical Hamiltonian mechanics (see, e.g., |l3)- 
Our approach clears up the presence of extra terms in the Law Hamiltonian and is equivalent to the method developed 
in the Law's paper [6|, at least, for quantum systems with quadratic Hamiltonian. 

Consider a massless scalar quantum field $(a;,i) in a nonstationary one-dimensional cavity formed by two mirrors. 
One of the mirrors is fixed at the point a; = 0, while the position of the other depends on time x = l{t). For 
definiteness, we impose the following boundary conditions at the mirrors: <I'(0,t) = <I>(/(t),i) = 0. 

The Hamiltonian for the field reads 




d<P{x,t)Y ^ (d<^{x,ty 
dx 



dx. (1) 



The field ^(x, t) and the canonical momentum H(a;, t) — d^{x, t)/dt are assumed to obey the equal time commutation 
relations [<I>(a;, t),H(a;', i)] — iS{x — x'), and hence the Heisenberg equations acquire the form 

<i)(x,t)-*[i/,$(a:,t)]=n(x,t), tl(x,t) =i[H, nix, t)]^^^^. (2) 

One can easily see that the Fourier coefficients 

Qn{t) = Wt^ / Hx,t)sm ( :^ 1 dx, (3) 




' n 



satisfy the standard commutation relations [Qn,Pn'] — iSnn' and [Qn,Qn'] ~ [Pn,Pn'] = 0. Introducing the "instan- 
taneous eigenfrequencies" uJn{t) = Trn/l{t), and the "instantaneous destruction and creation" operators a„, a|j 

which satisfy the canonical commutation relations [a„,a]j,] — Snn', [o.n,0'n'] = [O'niO.n'] ~ Oi one can write the 



expansions for the field and canonical momentum in the form 

'Trnx\ a„(t) + ajj(t) 




*(^'^) = Wt77T> sin 




(6) 

Since operators a„, ia\ and Qm Pn satisfy the same commutation relations, Eqs. lO can be considered as a canonical 
transformation. One can easily check that this transformation can be defined with a generating functional F[a^ <&] as 

^ dF[a,^ 5F[aM 

< = -'^^^ n(a:,t) = -^^^-^, (7) 

F[a, *] - -^ X! {^"WQ«[*' A+al- 2^2uJn{t)Qn['^, t]a„} • (8) 

n 

The order ambiguity in the last term of Eq. (jSJ can give rise only to an additive indefinite c- number constant, and 
thus is insignificant. 

For the stationary case {l{t) — const) the generating functional does not depend on time explicitly, so that the 
Hamiltonian remains unchanged under the canonical transformation (0),||SJ). In terms of new variables it acquires 
the form H — ^^ti;„ajja„. If the cavity is nonstationary, and hence the generating functional explicitly depends on 
time, the Hamiltonian changes, compare jlTj. The operators a„ satisfy the equations a„ ~ i[iJ',a„], where the new 
Hamiltonian is H' ^ H ~ dF/dt. 

For the case of a uniformly contracting cavity 

t <0, 
vt, Q<t<T, (9) 

t>T, 

where T — (If ~ li)/v is duration of contracting and u > is the velocity of the moving boundary, we obtain at 
< t < T 

H' = .mt). A - E H-r. - .Ual - af)} - ^ ^ (-1)"""'VW x ( (-"""' ;"»;"") + 4^1 . (10) 
^-^ y 47r J TT ^-^ 2m + n') rv — n \ 

The same result follows from the general Law formula (2.21) of Ref. fo] for the "effective" Hamiltonian which was 
constructed to reduce the equations of motion for operators Qn, Pn to the Heiscnberg form. This proves that both 
approaches are equivalent at least for the case when the Hamiltonian is quadratic in structure. 

The additional terms in H1U|) are responsible for annihilation and creation of pairs and scattering of particles from 
one mode to another in the contracting cavity. The corresponding "coupling constants" are proportional to the 
velocity of the moving boundary and are small in non-relativistic case. One can see that the terms of the Hamiltonian 
(|10|l responsible for interaction between modes are of the same order as those describing particle production in a 
separate mode. Therefore strong interaction between field modes is a built-in internal feature of DCE 18]. 

We see from Eq. (|10() that the Hamiltonian H' has a remarkable property. Its time dependence is determined 
exclusively by a c- number factor ir/lit). Hence, the operator A does not depend on time explicitly and is an integral 
of motion. Therefore, being quadratic in canonical variables a„, ajj, it admits diagonalization by a time independent 
linear canonical transformation (Bogolubov transformation) from a„, ajj to some new canonical variables 6„, 6jj. In 
terms of these new variables we would have 

A = E XJibn , H' = Y^i^^n/mPX- (11) 

n n 

Hereinafter it is convenient to introduce an auxiliary dimensionless time variable X = "" /n dt/l{t), which varies from 
— oo to +00. At < i < r, X = — (7r/w)ln(l — vt/li). In terms of this variable the Heisenberg equations of motion 
acquire the form 

da„/dx = i[A,a„], 0<t<T, (12) 

so that the operator A can be considered as the generator of translations in "time" x- The time moment t — 
corresponds to Xi = 0; whereas t = T to Xf = i''^/^') lii(l/p), where the dimensionless parameter p = If/U < 1 is the 
squeeze rate for the cavity. The value x = 00 corresponds to the moment t = li/v of the cavity collapse. 



III. DIAGONALIZATION OF THE HAMILTONIAN 



We will use the so-called Milne reference frame, see, e.g., Ref. [13, to diagonalize operator A. The corresponding 
spatial ^ and time r Milne coordinates are defined by 



t 



h 



rcosh^, x^rsinh^. 



(13) 



The point is that the cavity remains stationary in Milne reference frame. Indeed, the spacetime region Q < x < h ^ vt 
is mapped conformally by transformation (|13|) onto the strip < ^ < (1/2) In d, t > 0, where d = (1 + v)/{l — v) 
coincides with the Doppler factor for reflection from a mirror moving with velocity v. This means that the world line 
for the right boundary is given by ^ = (1/2) Ind = const. 
The wave equation in terms of Milne coordinates reads: 



1 a / a$ 

T dr \ dr 



1 92$ 



= 0. 



(14) 



This equation can be easily solved by separation of variables. Taking into account boundary conditions (which are 
stationary in the Milne reference frame) we obtain a complete set of modes 



which are normalized by the relation 



1 
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n=l,2,.. 



(15) 
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Note that the unusual sign in the LHS of this equation arises due to opposite directions of physical time t and Milne 
time r. In terms of the original variables x, t these modes acquire the form 



*„(a;,i) 



2\/wri 



_ v{t-x) \ >'"' _ / v{t + x) 



h 



h 



(16) 



and constitute a complete set of solutions for the wave equation in a uniformly contracting cavity orthonormalized by 
the usual Klein-Gordon scalar product 



m 

i I dx 

Hence, the general solution to the Heisenberg equations Q can be represented in the form 



d 
"9t ' 



<^{x, t)^J2 {^"*"(^' t) + bi^'Ux, t)} , 



(17) 



where operators 6„, 6jj satisfy the standard Bose-Einstein commutation relations and are independent of time since 
modes (I16f) satisfy the wave equation exactly. Specifically, this means that the Hamiltonian which governs behavior 
of the variables 6„ and bj^ is equal to zero with an accuracy to a c-number contribution. 

Let us now express operators 6„ in terms of operators a„, aj^. Due to completeness and orthonormality of the set 
((TE|l . we have: 



Ht) 



bn = i 



J ^n(.^,t)uix,t) - ^i^iilci>(x,t)| dx. 



(18) 



Using expansions © we reduce RHS of Eq. (|18|l to the form 

bn = J2 {^nn' {t)an'{t) + B„„, (i)a,V (i) | 



where 

m 



Now, using explicit expressions p6|) we can express time derivatives in Eq. (|19|) in terms of derivatives with respect 
to X and perform integration by parts. Introducing then a scahng variable y = x/l{t) wc finally obtain: 

Ann' (i) = e2™"^(*)/ '" -^ann' , B„„, (i) = e^^^xf*)/ '" '^/Sn™, , (20) 

where coefficients 

— 1 

- /dy(l-t;y)-2™VindgT«Wy^ (21) 
n J 

are already independent of time t, or auxiliary variable x- 
Thus, we have 

6„ = e2™"x(t)/in<i^^(^)^ (22) 

where operators 

init) = J2 {««"'«"'(*) + /9nn'4'(i)} . (23) 

n' 

are expressed in terms of operators a„, aj^ by a Bogolubov transformation with coefficients, which do not depend on 
time explicitly . Therefore the Hamiltonian which determines time dependence of the operators 6„ remains unchanged 
under transformation (|23|) and is given by Eq. p(J|) . This result, together with Eq. H22|) . is sufficient to ascertain the 
structure of the Hamiltonian A in terms of the operators 6„. 

Indeed, the Heisenberg equations of motion for operators 6„ are given by dbn/dx = «[A, 6„], compare (I12II . Taking 
into account the form (|22fl of time dependence of 6„, we have for the commutator [A, 6„]: 

[A, in] = -{2v/\nd)nbn. 

It clearly follows from pU|l and H23|l that A is quadratic in terms of operators &„ also. Hence, we unambiguously 
arrive to 

^ = |dE«^J^^«- (24) 

n 

Eqs. H24() . H23() and H21() constitute an exact solution for quantum field in a uniformly contracting cavity. This 
is one of the primary results of the present paper. Its physical meaning consists in existence of a notion of stable, 
non-interacting particles in a uniformly contracting one-dimensional cavity. Since pair production or mode-to-mode 
scattering are absent in this representation, the number of these particles in every mode is fixed. Their energies, 
however, depend on time due to collisions with the moving boundary, compare Eq. (|ll|l . It is worth noting that the 
exact energy spectrum of these particles turns out to be equidistant. 

It is important that the operators &„(i) for the case of a stationary cavity coincide with a„(i). Indeed, proceeding 
to the fimit u ^ in Eq. if^ . one can easily make certain that an,n'(0) = 6nn' , /9n,n'(0) = 0. We may conclude that 
operators 6„(t), bl^{t) defined by Eqs. (PS|l . (|^ represent a natural generahzation for particle destruction and creation 
operators valid for nonstationary cavities with uniformly moving boundaries. 

Let us now consider the Bogolubov transformation if^ with coefficients an,n'{v), /3„,„'(w) defined by Eq. if^ 
for v{t) being an arbitrary function of time. Since coefficients of the transformation (|23|l now depend on time, the 
operators 6„(i) are ruled by a new Hamiltonian H, so that the Heisenberg equations of motion acquire the form 

ibn = [H, bn]. Using Eqs. I|23(l and equations of motion for operators a„, we obtain for the time derivative of 6„(t): 



n = 2J ( "rrn'On' -|- /Jnn'fljj, -|- ann' O-n' + /?««'«„/ 
n' 

iann'an' + Prin' 0.;^, + «„„' «[i7', a„/] -|- /3„„/ i[i7', a,\,] j = ^ [ann'On' + ^nn'ol') +«[-ff'&n]- 



E 



6 

We know already the structure of H' in terms of 6„, see Eqs. (|10|l and (|24|l . For the difference H — H' = V we have: 

i[V, &„] = ^ iann'an' + Pnn'al,] . 
n' 

Using the well known unitarity relations for Bogolubov coefficients 

/ ^ \^nn" Pn'n" CiLn' n" Pnn" ) — U, / ^ \Otnn" ^n'n" Pnn" Pn'n" ) — ^nn' ^ V '^J 

n" n" 

which in our case (|21|l can be verified straightforwardly, we get the inverse transformation 

a« = X! y^n'rfin' " /?«'«&!/ j , (26) 



[V, bn]^vJ2 Unn'bn' + Bn^'bl) , (27) 



and hence, 



where 

•/vfin' — / ^ I , '-'n'n" Pn'n" i I J i-"?!™' — / ^ I , "n'n" Pn'n" , I • V^"7 

Substituting representations (|21|) for Bogolubov coefficients into Eqs. (j28|l we can perform summations over n" in H28|) 
using the well known relation 

OO (-. . oo 

^ ne""" = -- ^ 5'(x-2fc). 

n— — OO fc— — OO 



After this, all integrations become trivial and we obtain: 



In d , ^, 

+ 272(ln7)-l 

V 



TT (n' - n) (n' - n + i In d/27r) ' ' ' In^ d 

(29) 

/ 2")"'+" J ^^1 2-1ni(n' \n) I In d 



^nn' — ^" 



TT (n' + n){n' + n — ilii d/2Tr) 



where 7 = (1 — w^) ^/^. 

It follows from (f^ that .A*„, — —An'n, Bnn' — Bn'n- Using these relations and Eq. (|77|l one can easily find 

V = zvit) J2 {a.„' («(*))&' + I (Bnn'{v{t))Wn' ~ K^, («(i))6 J„,) | . (30) 

And finally 

^=I(^E-^t^" + V. (31) 

^ ^ n 

The representation (|31|l is very important because it clearly indicates that (i) the Hamiltonian for a field in a 
nonstationary cavity can be diagonalized if w = 0, and hence, we can use the notion of particles in a uniformly 
contracting cavity; and (ii) particle creation, as well as mode-to-mode scattering, takes place only when the boundary 
moves with acceleration. For our law of motion ^ these processes arise only in the initial i = and final t = T 
moments of time. The fact that radiation processes are absent during the time when one boundary of a cavity moves 
with a constant velocity was first noticed by Fulling and Davis in Ref. (,2Q|). 



IV. PARTICLE CREATION IN A UNIFORMLY CONTRACTING CAVITY 

We will now compute the number of particles created in a nonstationary cavity with the right mirror moving 
according to Eq. 1^1 • We will note first that approach based on the Hamiltonian H is not convenient for this purpose. 
This is because velocity of the mirror in our approximation changes instantaneously, and thus H contains ill defined 
products of 5 and ^-functions. Thereagainst, the Hamiltonian H' H1U|) is not that singular even in our approximation 
(it contains only ^-functions) and therefore operators an(t) are continuous functions of time. This point plays the 
central role for our method of solution. 

We start from the initial conditions 

a„(x)lx=-o = <", (32) 

where a™ are initial destruction operators which define the initial vacuum state by a™|0)i„ = 0. Due to continuity of 
operators a„(x) (a„(— 0) = a„(+0)), we have the initial data for the operators fon(x), 

^"(0) = Y. {«-'<" + /3n«'<"^} . (33) 

ri' 

According to Eq. (|22|1 the operator b^ depends on t (or x) by mean of phase rotation. SpecificaUy, dX t — T^ we have: 

&n(x/)=e-2-%(0), e = i^^. (34) 

Here the quantity Q measures the phase obtained by a particle in the principle mode (n = 1) during the squeezing 
period of the cavity. As it follows from Eq. (|34|l . the phase factor is completely determined by the fractional part of 
Q. For further convenience, we will denote the fractional part of Q by ^, 

^-™{^}- (») 

The range of variation of the quantity ?? is < i? < 1 (the symbol Frac denotes the fractional part of a number) . 

The operators anixf ~ 0)i which are the solutions for Eqs. H12() with the initial conditions H32() . can be now found 
from Eq. H26|l . Hence, again taking into account the continuity of operators a„, we finally have for a°"* = 
anixf + 0): 



out 



[an'J>n'{Xf) - f3n'nbi,{Xf)j ■ (36) 

n' 

By combining the three transformations H33fl . (|34|1 and (|36|) . we obtain: 
where the coefficients of the Bogolubov transformation (|37|1 are given by 

^ nn' — / ^ 1 ^ ^n"n^n"n' ^ Pn"n' Pn'^n r ; V*^^/ 

n" 

TA \ ^ J — 27rin"^ * /Q 27rm"^ * o 1 /''^Q^ 

^nn' — / ^ 1 ^ ^ri"nPri"n' ^ ^n"n'Pn"n f ■ V"^^ ) 

n" 

Consider the expression for Vnn' ■ Using Eq. ((2U, we get 



/. /. -1/1 \ 27rin'V In rf 

nn' '^ '^ 1/1 ..».- \ ' 



If f 1/1 \ ^~*'^ /ma 

- / dui / (kl'2 y^ — I ^^^ I X /g*''("J'i+"'a2)-27rm"iJ _ giTr(n'i/i+nj/2)+27rm"i?l /^q\ 

4 7 j J;^^ n" \l-vy2J I ' i ' 

After evaluation of the sum and double integral in H40|l (see Appendix for details) we obtain 

V , - ( n"+"' ^(C^-1)VW r -2™(d''-l)/(d-l) _ -27rin'id'-'>~l)/(d-l)\ /4-|^ 



The coefBcient [/„„' differs from y„„/ by the only change n' 
so that 



-n' everywhere, except for the common factor 



Un 



lT\n+n 



i{d^ 



27r(W-n')(n-n'di-'') 
FinaUy, the energy (mode) distribution of created particles is given by 



r -27rm(rf''-l)/(d-l) _ 2Tiin' (d^'^ -l)/(d-l)^ 



(42) 



Nnid,d) 



a: 



out '\ out 



'|0). 



EiK. 



{d-l)^n y-. 



°° n' sin^ 



(d 



-i^{n + n'd'-^) 



{nd'^ + n')^{n + n'd^-^y 



(43) 



One can see that the number of particles H43|) created in each mode depends on two parameters. One of them is 
the Doppler factor d, which can be expressed either in terms of velocity of the moving boundary, d = (1 + u)/(l — v), 

or in terms of the corresponding Lorcntz factor 7 = (1 — w^)^^/^, rf = ( 7 + vT^^l) ■ The second parameter 1} is 

determined by d and a dimensionless squeeze rate of the cavity p — If /U which, in turn, can be expressed in terms of 
the duration of squeezing T as p = 1 — vT/li. Note that acceleration of the moving boundary. 



W 



-vS{t)+vdit-T), 



is determined by the same set of parameters. It is natural, since particles can be created by an accelerating boundary 
only. 

An interesting feature of the obtained solution is that the number of created particles is a periodic function of the 
parameter log p with the period log d. This can be traced already from the expansions (|38|) , (|39|) , and definition l|35|) 
of the parameter i?. Surprisingly, the number of created particles completely vanishes if i9 = 0, see Eq. (|43|l . It is 
worth noting that Nn tends to zero also if t? ^ 1, so that the energy spectrum of created particles is a continuous 
function of parameter p. According to Eq. (|35() . 1? turns into zero, and hence the number of created particles vanishes. 
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FIG. 1: A graphical illustration for null values of the number of created particles at a) p : 
represent world lines of a light ray, emitted by the moving boundary at t = 0. 



1/d, b) p = 1/d'. The dashed lines 



Up— 1/d^, where fc = 0, 1, 2, . . .. The value fc = corresponds to a stationary case, while the others can be explained 
by destructive quantum interference between events of particle creation at the moments of acceleration {t — 0) and 
deceleration {t = T) of the boundary. As it is illustrated by Fig. ^ the values of the squeeze rate p = 1/d^ exactly 
correspond to the cases when a light ray, emitted by the right boundary at t = 0, returns to it at i = T after 
fc — 1 successive reflections and therefore can be absorbed at the moment of deceleration. Analogous arguments were 
employed by the authors of Ref. (i2(j|) in discussion of the effect of vanishing of the renormalized energy flux in a 
uniformly expanding cavity, though only the case fc = 1 was considered in that work. It is worth mentioning that at 
p — 1/d^ every level of the initial energy spectrum of the fleld inside the cavity converts in time interval T to the 
value w^ = nn/lip = {nn/li)d'' coinciding with the frequency acquired by a particle after fc — 1 successive reflections 
from the moving boundary, if this particle was emitted by the boundary moving with velocity w at i = 0. 

The average number of particles, created in the principle {n = 1) and in the first excited (n = 2) modes as functions 
of 1} are shown in Fig. |5] for several values of the moving boundary 7-factor. The same plots after the appropriate 
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FIG. 2: The average number of particles created in a) the principle mode (n = 1), and b) the first excited mode (n = 2) as a 
function of i9 for different values of the moving boundary 7-factor. 



scaling represent the first period of 7V„ dependence on log(l/p). 

One can see from Fig.|21that at sufficiently large 7 plots for both N\ and iV2 have a plateau at not too small ?9, and 
a strongly pronounced maxima (one for n = 1 and two for n = 2) at ■(? close to unity. It is seen also that the altitudes 
of the plateaus are growing with 7 remaining however rather small. To estimate the maximum possible altitude of 
a plateau we will note that in the limiting case 7 — > 00 (d — > 00) at a fixed p the instantaneous approximation 
considered in Ref. 13\ becomes applicable. Indeed, one can easily make certain that Eq. H43|l reduces to Eq. (6) of 
Ref. 13] under condition d 3> n(l — p)/ p ^. For our estimation we need only high-energy tail of energy distribution 
in instantaneous approximation (see Eq. (10) of Ref. [13I) where we should put p — since not too small values of d 
at d 3> 1 correspond to small p. Thus, for maximum values of 7V„ in the plateau region we have 



ivip') 



1 



2n'^n 



[ln(27m) + C - 1] 



(44) 



where C = 0.5572 is the Euler constant. For n = 1, 2 we have N^^' « 7.17 • 10"^ and N^^' « 5.34 • 10"^ respectively. 
So, in the plateau region the average number of created particles is small even for ultrarelativistic velocities of the 
moving boundary. 

To estimate positions of the maxima and the peak values of Nn we should evaluate the sum in Eq. H43|l under the 
conditions d ^ 1, 1 — •& -^ 1. It is clear that the major contribution to the sum is given by n'^ff ^ nd"^ . Thus we 
have: 



00 , -2 



Nn 



E 



n' sin^ [7r(W-i - n'd-^)] 



7r2 ^^ (n + n'd-'')2[n' + W-i] 



4- sin^ (nnd^-^) V — 



'^ (n' + nd''-i)2 



Since nd^ 3> 1, we can change summation over n' by integration. Finally, after evaluation of the integral we obtain 
in the leading logarithmic approximation: 



Nn « ^^ sin^ (nnd'^-^) \n{nd). 

We see from Eq. ()45|l that positions of the maxima are determined by the relations 

7rnd''-i = ^(j + 1/2), J = 0, 1, ... (n - 1). 



(45) 



(46) 



^ Note, that this condition of applicabiUty of instantaneous approximation was obtained in Ref. Il3l on the basis of merely quahtative 
arguments. 
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This is an approximate condition but it becomes exact in the ultrarelativistic hmit. For example, at 7 = 2, 10, 10^, 10'^ 
for the single maximum position at n = 1 we have respectively dmax = 0-74,0.88,0.94,0.95 from Eq. (|46|) . and 
i^max = 0.85, 0.90, 0.94, 0.96 as the result of computation according to exact formula H43|l . 
The condition (|46|l can be reformulated in terms of wavelengths as 

{2j + l)Xfj2^2k/d'^+\ (47) 

where A4 — '2lf jn is the wavelength corresponding to the n-th harmonic in the final state. The physical meaning of 
this condition can be clarified by the following qualitative consideration. 

We have shown, see Eqs. (I31|l . l|30() . that particles are created inside the cavity only at the moments of start [t ~ 0) 
and stopping (t = T) of its boundaries. Physically, it is clear that in the framework of the applied approximation 
(instantaneous acceleration) the created particles constitute an extremely narrow wave packet (cluster) which prop- 
agates then inside the cavity without spreading successively colliding with its boundaries, see Fig. 13 Such picture 
agrees with the results of calculation of the energy momentum tensor performed by Fulling and Davies in Ref. (|2^). 

Let us describe the clusters of particles created at the moments t = and t — T hy classical fields <i>i(a;,i) and 
<i>2(a;,i) respectively. The function <i>2(a;,r) describes the second cluster produced by the boundary stopping at the 
moment t = T. It is located exactly at the point x = lf. The first cluster at the same moment of time is located 
at some distance A left to the first one. The form of the cluster $1 is absolutely the same and differs from $2 
only by sign since the boundary accelerations at two moments of particle creation were of opposite sign. Thus, 
^i{x, T) ~ —$2(2; — A, T). Expanding the fields <i>i in terms of stationary owi-modes 

Vr (^, = -^ sin f ^x) e-*"/"*, UJfn = {^/lf)n, 

for the total field $(x, T) = <^i{x, T) + $2(2-, T) we obtain 

Thus, we see that contribution of the n-th mode to the field $(a;, T) is absent if nA/2lf — j, j — 0, 1, 2, ... (destructive 
interference of two clusters) , or it is maximal if 

nA/2lf=j + l/2, j = 0,1,2,... (48) 

(constructive interference). In view of the restriction A < /j, we have 2j + 1 < n. Since particle creation and mode 
interaction are absent in a stationary cavity, these conditions hold at arbitrary moment of time t > T. Now, the 
question is whether it is possible to match Eq. H48() with conditions H47|l (as well as the corresponding conditions for 
minima) . 

Assume that the first cluster undergoes k collisions with the moving boundary during the period of contraction T. It 
can be easily found that the last collision takes place at the point Xfc = k/d'' at the moment of time ife = {li/v){l—l/d''). 
During the period of time T — tk the cluster either collides with the boundary at rest, if T — i^ > Xk (Fig.|3^), or not, 
if T — <fe < Xk (Fig. Ob)- The quantity A for these cases is given by the following relations respectively: 

A^lf±[xk-{T-tk)]. 

Then Eq. (|48ll is equivalent to two relations 

T -tk> Xk 

(49) 
T -tk <Xk 

where j < {n — l)/2. In the ultrarelativistic limit 1 — v ^ 1 {d ~ 2/(1 — v)) the conditions H49() take on form of 
Eq. H47|l with < j < n — 1. Our analysis shows that physically the conditions H47|l mean that maxima of Nn appear 
in the cases when the first cluster is located at a distance equal to an odd number of half-waves X(^/2 from the point 
of creation of the second cluster. It is worth emphasizing that kinematics of our problem permits implementation of 
this condition only in the ultrarelativistic case. Therefore the less is the speed of the boundary u, the less pronounced 
are the maxima and N^ 7^ in minima, compare Fig. |21 

The values of Nn in the points of maxima are given by Eq. (|45|l with logarithmic accuracy, and thus cannot reproduce 
numerical results just as well as their positions. However, Eq. (|45|) reproduces the main features of spectra nearby the 





( 2lf f 1 + v . 1 


h l~v 


1 n r 2« ^2 


d^ V 


2lf f 1 1-v 

\ 1 \ 1 Tt 




^ n V 2 V 
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a) 



.t A 




K 


1 


X 




FIG. 3: World lines of cavity boundaries (solid lines) and of the cluster of particles (dashed line) created at f = for a) 
T — tk > Xk, and b) T — tk < xt- The number of collisions with the moving boundary k = 1. 

maxima qualitatively correctly. It is important that the height of the maxima grows as log d, or equivalently, as log 7 
in the ultrarelativistic limit due to strong interactions between the modes. Clearly, the effect of DCE amplification 
arises due to constructive interference of events of particle creation at t = and t = T . For example, if 7 = 10"^, then 
for certain values of p (or T) in the mean more than one particle in the principle mode (or a little bit less than one 
particle in the first excited mode) can be created. 



DISCUSSION 



An exact solution to the problem of dynamical Casimir effect in a one-dimensional uniformly contracting cavity 
was found. This became possible due to existence of stable non-interacting particle states in such a cavity. As a 
consequence, we succeeded in diagonalization of the Hamiltonian for a many-mode quantum field in a cavity with 
uniformly moving boundaries (a generalization to the case with both boundaries moving uniformly is straightforward, 
but rather cumbersome). The developed method opens us possibilities to solve exactly many other one-dimensional 
DCE problems for cavities with boundaries moving along polygon-like world lines, including the interesting case of 
vibrating boundaries. A new representation for the field Hamiltonian was proposed in one-dimensional problem. It 
was demonstrated explicitly that particle creation as well, as mode-to-mode interaction is absent in a nonstationary 
cavity if its boundaries are moving uniformly. However, the same approach can not be applied to physically more 
realistic three-dimensional problem, although a complete set of classical solutions generalizing Eq. (|15|l is known for 
this case as well l9| . This is due to breakdown of stability of particle states and presence of interactions between them 
in three-dimensional problem. 

It was shown that the number of particles, created in each mode, depends periodically on the logarithm of the 
squeezing rate of the cavity p, which is equal to the ratio of the final and initial sizes of the cavity and is related to 
duration of contracting in a rather simple way. This dependence is not monotonous. For the given velocity of the 
moving boundary, there exists an infinite sequence of values pk for which particles are not produced. On the other 
hand, it is possible to specify an infinite set of combinations of parameters, which correspond to an optimal regime 
of particle production. Theoretically, at sufficiently large Lorentz factor 7 for the moving boundary, the number of 
created particles can be done arbitrary large. It is worth noting, however, that the mentioned values of Lorentz factors 
are indeed very large. For example, the value 7 ^ 10"^ for production of a single particle is necessary. We should also 
mention that the number of created particles in the case of a cavity with ultrarelativistic boundaries depends on 7 
rather slowly (proportional to log 7). 

The non-monotonous character of dependence of the number of created particles on p is also present in the frame- 
work of an oversimplified "one- mode" model. However, correct estimations, especially near the maxima of particle 
production rate, and even correct qualitative understanding require taking into account mode-to-mode interactions. 
This is true in principle at consideration of any non-perturbative feature of DCE. Therefore, it is very important to 
consider exactly soluble models. It should be noticed, that the approximation of uniform motion of the boundaries, 
separated by moments of instant infinite accelerations, does not work well for high-frequency modes. For example, 
the energy spectrum (|43|l has an inadmissibly slowly decreasing tail (7V„ oc 1/n) leading to an infinite total number of 
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created particles. These problems can be resolved by smoothing instant infinite accelerations, but the corresponding 
problem ceases be exactly soluble. 
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APPENDIX A: EVALUATION OF THE INTEGRAL IN EQ. l|40|l 

The coefficient Vnn' (|40|) is represented by a difference of two integrals. After the change of integration variables 
2/1 ^^ 2/2 in the second of them, T^„' reduces to the form: 

2Trin"Z{yi,y2) 



ly nn , , i^/„, ,,_!_„', ,„^T x ^ e' 



where 



dyi / d2/2e^"^"^^+"'"=^Ini ^ , (Al) 



1 -1 ""^1 



-^ 1 - vyi 



Z{y,,y2) = log, d-^- ^ . (A2) 

1 — vy2 / 



For calculation of the sum in Eq. (|A1|I we will use formula 1.441.1 from |23| generalized for arbitrary values of Z 



Im ^ - — j^— = 7r[--z)+7r ^ ke{Z - k)e{k + I - Z) . (A3) 



n">l 



Obviously, the first term in the RHS of ljA3|l gives zero contribution to the integral (|A1|) . Taking into account that 
— 1 — "(9 < Z[yi, j/2) < 1 — ^ inside the integration domain, and < ^ < 1 by definition, we obtain 



1 1 

K„' = -^^^ fdy, I dy2e'<''y^+^'y-'^{e{-z)e{z + i) 



(A4) 



-1 -1 

+2e{-z -i)e{z + 2)}. 

It is clear that the combination of ^-functions 

e{z) + e{-z)e{z + 1) + e{-z - i)e{z + 2) 

is equivalent to 1 in the domain of integration, and thus 

^^2y g»7r(nyi+«'a2) ^Q(^z) + e{-z)e{z + 1) + e{-z - i)e{z + 2)} = 0. 

Therefore, we have: 

K™' = '-^^^ jj d'y e-("^i+"'^^) {e{Z) -e{-Z- 1)} . (A5) 

The integral IjASp can be easily reduced to the form 

1 1 c{Vl) 

i-Ky/nn' 



y„„, = l^i^ S / '^yi y dy^ -JdyiJ dy, } e"("^^+"'^^) , (A6) 

-1 6(yi) a -1 
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where a = 1- 2{d^ - l)/{d - 1), 6(yi) = (1 - d-^)/v + d'^yi, and c{yi) = (1 - d^-'^)/v + d^-^yi. After simple 
calculations we now obtain the expression H41|) for Vn.n'- 
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